Electroweak Bremsstrahlung in Dark Matter Annihilation 
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A conservative upper bound on the total dark matter (DM) annihilation rate can be obtained 
by constraining the appearance rate of the annihilation products which are hardest to detect. The 
production of neutrinos, via the process XX ~' '^^y has thus been used to set a strong general 
bound on the dark matter annihilation rate. However, Standard Model radiative corrections to 
this process will inevitably produce photons which may be easier to detect. We present an explicit 
calculation of the branching ratios for the electroweak bremsstrahlung processes xX ~^ vvZ and 
XX ~* veW . These modes inevitably lead to electromagnetic showers and further constraints on the 
DM annihilation cross-section. In addition to annihilation, our calculations are also applicable to 
the case of dark matter decay. 

PACS numbers: 95.35.+d, 95.85.Ry 



I. INTRODUCTION 

The identity of the dark matter (DM) is one of the 
great unresolved questions in particle physics and cos- 
mology P, 0, Q . An important method of probing DM 
properties is via indirect detection, whereby we look for 
the appearance of particles produced via DM annihila- 
tion or decay. We can search for such a signal ema- 
nating from the dark matter concentration in our own 
galaxy, other galaxies (satellite, dwarf, or clustered), or 
for a isotropic flux from the dark matter distributed 
throughout the Universe [J] . Investigated signals include 
positrons, gamma-rays, x-rays, and even microwaves (the 
"WMAP haze" 0). 

If we make the reasonable assumption that DM decay 
or annihilation products must be Standard Model (SM) 
particles (i.e. we assume the dark matter is the lightest 
stable particle in the beyond-SM sector) then it is pos- 
sible to set a conservative upper bound on the total DM 
annihilation rate by looking for the annihilation products 
which are the hardest to detect, namely, neutrinos [6|. 
All other possible final states would lead to the produc- 
tion of gamma rays, for which more stringent bounds 
apply. For example, quarks and gluons hadronize, pro- 
ducing pions and thus photons via 7r° — > 77; the decays 
of r^, ^ and also produce ifi . Charged parti- 
cles produce photons via electromagnetic radiative cor- 
rections 0, while energy loss processes for also 
produce photons @. By calculating the cosmic diffuse 
neutrino flux produced via the DM annihilation process 
XX ill sill halos throughout the Universe, a strong 

and general bound on the DM total annihilation cross 
section has recently been derived @ . The corresponding 
signal from our own galaxy can be used to set a compara- 
ble limit (and improves upon the cosmic bound in some 
mass ranges) [l^l while the technique has been extended 
to low (MeV) masses in Ref. [Hj . Analogous bounds have 
been derived for the DM decay rate 12]. 

The general upper bound on the total DM annihilation 
cross section defined via the limit on xX ~^ is surpris- 
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FIG. 1: The lowest order tree level process xX ~^ 
(left) is accompanied by electroweak bremsstrahlung pro- 
cesses (right). 



ingly strong. (See Ref. [Tjl for a comparison between 
photon-based and neutrino-based limits.) However, a 
scenario in which neutrinos alone are produced in the 
final state is technically impossible. Even leaving aside 
the theoretical issue that a direct coupling of DM to only 
neutrinos violates the S'?7(2)-invariance of the weak inter- 
action, electroweak radiative corrections imply indirect 
couplings to states other than neutrinos. For example, 
for energies above My/^z-, electroweak bremsstrahlung of 
W or Z bosons can occur at sizeable rates [13, [ll] , see 
Fig. [TJ The hadronic decays of these gauge bosons pro- 
duce neutral pions, which decay to gamma rays. Even 
for energies below Mvi/, processes involving virtual elec- 
troweak gauge bosons will lead to particles with electro- 
magnetic interactions, though the rate for such processes 
is suppressed at low energy. 

Kachelriess and Serpico have estimated the constraints 
on the cross section for xX ~^ '^'^ (and hence on total DM 
annihilation cross section) by considering gamma ray s 
produced via the accompanying process xX ~^ vvZ [14| . 
We present here an explicit calculation of the branch- 
ing ratios for the electroweak bremsstrahlung processes 
XX ~* vvZ and xX ~* veW . 

Note that the expected magnitude of the DM total an- 
nihilation cross section varies enormously between spe- 
cific models. For LSP (lightest supersymmetric parti- 
cle) DM, s-wave annihilation to fermions is helicity sup- 
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FIG. 2: Feynman diagram for the W-strahlung process XX 



pressed and thus the lowest order annihilation rate can be 
quite small. On the other hand, there is no such suppres- 
sion for Kaluza-Klein DM. Please refer to the Appendix 
for a detailed discussion of these issues. Our analysis is 
not specific to any particular dark matter candidate (nor 
do we assume the DM is a thermal relic). 



II. W-STRAHLUNG 

We shall first consider DM annihilation to four body 
final states via W-strahlung, an example of which is the 
XX ~* e^e~vv process shown in Fig. [21 and calculate the 
ratio of the cross section for this process, to that for the 
lowest order tree level process xX ~^ For simplicity, 
we will assume that the coupling between the DM-current 
and the neutrino-current is mediated by a scalar boson 
"B" . Given this scalar coupling, the terms in the matrix 
element involving the initial state (x) particles will fac- 
torize from the full matrix element. It is thus useful to 
consider the matrix element for the decay of the virtual 
B*. 

We will first discuss the phase space calculation, and 
begin by noting that it is useful to treat the process 
shown in Fig. [2] as a sequence of three 1 — > 2 particle 
decays of a virtual particle, viz., B*[P) —f y*{q) + v{ki), 
followed by v*{q) W*+{Q) + e-{k2), followed by 
W*+{Q) v{l2) + e+{h). The four-body final state 
Lorentz-Invariant Phase Space (LIPS) is given by 

LIPS^'^^ = {27r)^ J dk2 J dki J dh J dh 

X S^(P-h-l2-k^~k2), (1) 

where dk = {2Tr)^^ d'^k/2ko, etc. We integrate over the 
momenta of the virtual particles by inserting 

dUd^Q6\q-k2-Q)S\Q-h-l2), (2) 

which ensures momentum conservation for the virtual 
processes. The phase space then factorizes into a product 
of three separate two-body phase space factors convolved 
over the two virtual particle momenta: 

LIPS^^^= I ^1^1 dLIPS^^\P\q\kl) 

X J dLIPS^^'^{q'^,Q^,ki) J dLIPS^^\Q^,llll). (3) 



Each two-body differential phase space factor is easily 
evaluated in the respective two-body center of momen- 
tum (CoM) frame, using the expression 



<iL/PS<«(.,„,=).i^fe3(^,, (4) 



where dH. is the CoM solid angle, and X{x,y,z) is the 
triangle function given by 



\{x, y, z) = + y"^ + — 2xy — 2xz — 2yz . 



(5) 



The phase space may be written in a form useful for our 
calculation. 
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^) JdLIPS^'^iQ'AO). (6) 



Here and throughout, we neglect the masses of the lep- 
tons, hence fcj = = = = 0. 

We next calculate the matrix element for B* 
De e~W~^* — > De e~h'e e"*", which is given by 



Mw = gs- 



u{k2)j' 



.1-7' 



1 _ 

w(^2)7m — o — "(^i' 



■— Jt'(^l) 

q^ — mf, 
-1 



(7) 



where the (non-standard) Biyv and (standard) Wve cou- 
plings are gs and g/V2 = e/\/2 sin^u,, respectively. We 
have expressed the matrix element in Eq. ^ in Feynman 
t'Hooft gauge, but note that our result is gauge invariant. 
After squaring, summing over spins, and performing the 
integration over the momenta li and I2, we obtain 



w\ 



Ss {Q ■ k2){Q ■ q) 



4s (fe -g) 
1 



spnis 

f {Q ■ k2){Q ■ P) ~ 4Q\^ {k2-P) 

7 4 

9b9_ 

3-26 (27r) g4 [(Q2 _ ^2^)2 + ^mwTw?] ' ^ ^ 

where s = is the center of mass energy squared, and 
we have made use of the identity 

[ dLIPS'^^^Q^, 0,0) I'^l^ - [2Q"Q'3 + Q2g"/3] . 
J 9o7r 

. 

We now multiply this partial result by the remaining 
part of the phase space in Eq. ([S]), and perform the 
integrations over dfip, d^lq, and q^, to obtain the rate 
r(— > i>ei'ee^e~). We wish to compare this rate with that 
for the B* v^Ve, for which the lowest-order tree- level 



expression is / dLIPSJ^ 
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gls/iAn). The 
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resulting expression for the ratio of rates for these two 
processes is 



r(^ Vee'W*+ VeVee^e-) 
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■^w (10) 



r(^ u.D,) 32 29 (27r)4 

" V 1 + Uy^ - 9y2 -9y- (6y^ + ISy"^) ln(?/) 
X j ay 

'0 



expressed in terms of dimensionless scahng variables y ■ 
Q^/s and xw = s/m^. 
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FIG. 3: Integrand of Eg . [TOl versus y 
order of the curves, the values of xw - 
1.1, 1.5, 2.0, and 5.0. 



Q^/s; in ascending 
s/m^ are 0.9, 1.0, 



In Fig. [3] we plot the integrand of Eq. (fTO|l versus 
y = Q^/s, for the values xw = s/m^Y — 0.9, 1.0. 1.1, 
1.5, 2.0, and 5.0. The figure reveals that, for xw ^ 1-5, 
i.e, for nix ^ ^ O.Gmw, the rate r(— > I'ei^ee'^e^) 

is dominated by the on-shell M^-resonance. (Apparently, 
the q^^ pole in Eq. ([5]) is effectively negated by the van- 
ishing of massless three-body phase space for 1/ e+e^ at 
= 0.) Thus, we are justified in using the narrow width 
approximation (NWA) for the W^-propagator 



(5(Q2 _ m2). (11) 

AIT ) \ ) 



(Q2 _ ^2) _^ ^MTf 



With this approximation, the cross section factorizes into 
the on-shell production and subsequent decay of the W 
boson, and the contribution from virtual (off-shell) W 
bosons is neglected. Using the NWA approximation, the 
integral in Eq. pUj) is easily evaluated. The resulting 
ratio of widths becomes 
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(12) 
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Here we have dropped the flavor subscript on the 
(anti)neutrinos since this ratio remains the same when 
flavors are summed, and multiplied by a prefactor of 2 x 9, 
which we now explain. The "2" comes from adding the 
VF~-strahlung channel to the channel. Note that the 
two amplitudes do not interfere since the charges of the 
produced M^'s, and therefore of the pairs they produce 
with invariant mass M mw, are distinguishable. The 
"9" comes from summing over all decay channels avail- 
able to the decaying W. We have three leptonic channels, 
and two quark flavor channels, the latter multiplied by 
three for color channels. 

We may also evaluate the W-width. At the level of our 
calculation, we have for this quantity 



Tw = (9) X ^ mw ■ 



(13) 



The "9" here is the same final state count that appeared 
in Eq. p^ . Inputting this width into Eq. and using 
— 47rQ!/sin^ 9^, we arrive at our final expression. 
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V 48 / 



(14) 
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We note that this expression may also be obtained by di- 
rectly computing the production of real (on-shcU) gauge 
bosons. However, in that case one must choose unitary 
gauge, where all degrees of freedom are physical, in order 
to reproduce Eq. p4)) . 

For our numerical work, we will take sin^ 9^ = 0.231 
[sm^"^ 9w — 4.33), and a — 1/128 as appropriate for 
physics at the electroweak scale. This latter choice is es- 
pecially appropriate in light of the accuracy of the NWA; 
q^ of the virtual neutrino will itself have a value near the 
threshold for on-shell W production, i.e., at ^ m^. 



III. Z-STRAHLUNG 

The cross section for the Z-strahlung process, xx ^ 
vvZ* All, may be calculated similarly to that for W- 
strahlung. For a scalar coupling as assumed in Fig. [21 
there is no interference between diagrams in which the 
Z is radiated by the v and v. Thus the cross section 
for the Z channel is simply obtained from that for the 
W channel given in Eq. by dividing by a factor of 



2 cos^ 9 



w 



1.54, and replacing xw with xz = s/m 



2 



Rz{xz) 



1 



2 cos^ 9w 



Rw{xz) 



(15) 



IV. DISCUSSION 



Figure [3] shows the ratios Rw and Rz as functions 
of my. We choose to plot rates versus to^ rather than 



4 



xw and xz to make the presentation more physical, and 
to better illustrate the difference between the W and Z 
rates. To convert from the scaling variables to m^, we 
have used the expressions xg = s/mQ « 4:{m'^/mQ), 
G = W, Z , appropriate for non-relativistic dark matter. 
The curve for Rz may be directly compared to the points 
from [14], which we show. Qualitative but not quantita- 
tive agreement is evident. 
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FIG. 4: The ratios of rates for W- and Z-strahlung to 
production are plotted versus the dark matter mass m^. Evi- 
dent is the the dominance of the leading linear term in Eq. [121 
above xq ~ 10. Extrapolations beyond ~ TeV hold some 
uncertainty due to multiple bremsstrahlung, to possible yet 
unknown new physics, and eventually to re-summation of in- 
frared near- sing ularities. Shown for comparison are the Rz 
points from [14| . 

Let us discuss some general features of Rg, G = W,Z. 
At large s 3> itiq , we expect a leading term linear in the 
dimensionless variable xg = s/ttiq. The factor of s arises 
from the ratio of 3-body to 2-body phase space, while the 
numerator is provided by the only other dimensionful 
quantity in the process. From Fig. |4l we see that the 
leading linear term indeed dominates above xg 10, 
which corresponds to > 1.5 m^. Thus, we may write 
a very simple expression for the width ratio at xg ;J 10. 
It is 




for xw>10, (16) 



and likewise times (2 cos'^^h')^^ for Rz- It is unsurpris- 
ing that the inequality s > 10 m^. 4^ > 1.5 mc has 
appeared twice, earlier to put the W oi Z on-shell, and 
here to impose the dominance of the leading term in the 
expression for W- or Z-strahlung process. 

In the very large s (or equivalently, the very large m^) 
limit, the branching ratio for multi W/Z production will 
become sizeable. We may estimate the onset of double- 
W/Z production. The general formula for n-body mass- 
less phase space is 

^^^^ " " 8^ V 16^J (n-l)!(n-2)! ' ^^^^ 



Neglecting combinatoric factors, the perturbative ex- 
pansion parameter for additional W/Z bosons is then 
(-^)^ ~ am'^/mQ. Thus, perturbation theory be- 
comes unreliable and multiple W/Z production occurs 
when Xg ~ 2T:/a, or ^ a~^^^ tug ~ TeV. Resumma- 
tions in the very large s regime, involving ordered ln^(a:G) 
terms from emission of (nearly) massless or coUinear 
W/Z's, are discussed in [l^. Co-emission of a hard pho- 
ton will also occur, at a rate comparable to double W/Z 
emission [l6l |. 

Finally, we note that our results are easily applied to 
the DM decay process x ~^ i'^? where x is now a boson, 
with Xg = s/rriQ —f m'^/mQ. Indeed, similar expressions 
will hold for any neutrino production mechanism in which 
the invariant mass of the Dv exceeds mw- 



V. CONCLUSIONS 

The decay of W and Z bosons produced via elec- 
troweak bremsstrahlung will lead to neutral pions and 
thus photons. One may constrain this DM annihilation 
signal by considering its contribution to the Galactic or 
extragalactic diffu se g amma ray background. This was 
considered in Ref. [ij] , where it was shown that the con- 
tribution of the process a-{xx ^ vvZ^ to the galactic 
gamma ray background imposed limits on the lowest or- 
der process, (j{xx ~^ comparable to those obtained 
directly with neutrinos. The branching ratio expressions 
we have derived differ quantitatively, through not qual- 
itatively, from the cross-section estimates in Ref. _l4 ] . 
Thus, our results lead to similar bounds. Future Galac- 
tic gamma ray observations, such as those to be made 
by GLAST, have the potential to somewhat reduce the 
diffuse backgrounds through better point source identifi- 
cation, and to measure the background more precisely. In 
turn, this will strengthen the electromagnetic constraint 
on DM annihilation, and increase the utility of the quan- 
titative results we have presented herein. 
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VI. APPENDIX: INITIAL (xx) OR FINAL (uu) 
STATE MAJORANA FERMIONS 

Although our goal in this paper has been to present 
radiative corrections to the no-neutrino tree-level pro- 
cess in a manner as model-independent as possible, it 
is nevertheless interesting to ask what constraints would 
arise if the particles of either the intial xX state or the 
final i>u state (or both) are Major ana particles. It is 
quite possible that neutrinos are Majorana particles. It 
is also possible that the DM is Majorana. Although the 
LSP in supersymmetric extensions of the SM is not a 
neutrinos-only DM candidate, it provides a popular ex- 
ample of Majorana DM. SUSY examples of Majorana 
fcrmions include the ncutralino and the photino. (On 
the other hand, Kaluza-Klein DM, with mass reflecting 
the length scale of extra dimensions, is a popular example 
of non-Majorana DM. Typically the LKP is the bosonic 
recurrence of the photon. For a no-neutrinos model, the 
DM would be different again, and currently unknown.) 

Two identical fermions comprise a Majorana pair. A 
fermion pair can have total spin S in the symmetric state 
S' = 1 or in the antisymmetric state S ~ 0. The par- 
ity of the two-fermion state is P = {—)^~^^, where L is 
the orbital angular momentum of the pair. This parity 
formula holds for both Dirac and Majorana pairs. The 
negative intrinsic parity of the pair, independent of the 
orbital parity (— )^, is the same for Dirac and Majorana 
pairs for different reasons. In the Dirac case, the u and 

V spinors (equivalently, the positive and negative energy 
states) are independent and have opposite parity corre- 
sponding to the ±1 eigenvalues of the parity operator 7^. 
Reinterpretating the two spinor types, or positive and 
negative energy states, as particle and antiparticle, then 
leads directly to opposite intrinsic parity for the particle- 
antiparticle pair. In the Majorana case, the fermion has 
intrinsic parity ±i, and so the two-particle state has in- 
trinsic parity (±i)^ = —1. 

On general grounds, the L*'' partial wave contribution 
to the annihilation rate is suppressed as w^^, where v is 
the relative velocity between the heavy, non-relativistic 
XX pair. The virial velocity in our Galactic halo is only 

V ^ 300 km/s ^ lO^'^c, so even for L = 1 the supression 
is considerable. Thus only the L = partial wave gives 
an unsuppressed annihilation rate in today's Universe. 
The L > 2 states are too suppressed to contribute to 
observable rates. 

A Majorana pair is even under charge- conjugation 
(particle- antiparticle exchange), and so from the general 
relation C = {—)^^^ — +1 one infers that L and S must 
be either both even, or both odd for the pair. The ori- 
gin of the C = (— )^^'^ = rule is as follows: Under 
particle-antiparticle exchange, the spatial wave function 
contributes (— )^, and the spin wave hmction contributes 
(-1-1) if in the symmetric triplet S = 1 state, and (—1) if 
in the antisymmetric S = singlet state, i.e., (— )'^+^. In 
addition, there is an overall (—1) from anticommutation 
of the two particle-creation operators b^d^ for the Dirac 



case, and b^b^ for the Majorana case. 

Consider the L < 2 states. In spectroscopic nota- 
tion (^'^+^)Lj and spin-parity notation (J^^), the vector 
'^Si (1 ), C-odd axial vector ^Pi (1^ ), and assorted 
'^Dj (J~~) states are all C-odd and therefore disallowed. 
The pseudoscalar ^Sq (0 '"), scalar ^Po (O^"*")? axial vec- 
tor ''Pi (1^^), C-evcn tensor ^P2 (2++), and pseudoten- 
sor (2 ^) are all C-evcn and therefore allowed. In 
particular, the sole L = state, with no v"^^ suppression, 
is the pseudoscalar ^Sq (0 

Incidentally, at threshold, defined by s = {2M^)'^ or 
V = \J^ — 4-^1 /s = 0, the orbital angular momen- 
tum L is necessarily zero. With two identical Majorana 
fermions, the two-particle wave function must be anti- 
symmetric under particle interchange. Since L = at 
threshold, the xX spatial wave function is even, and the 
wave function must be antisymmetrized in its spin. The 
antisymmetric spin wave function is the 5 = state. 
Thus, the only contributing partial wave at threshold is 
the state. We have just seen that this is also the 
only state with no v"^^ suppression, so one may expect 
an unsuppressed Majorana annihilation rate at threshold 
if and only if there is a ^Sq partial wave. 

One may also invoke CP invariance to note that the 

spin S of initial and final two-fermion states, Dirac or 
Majorana, are the same. This follows simply from CP = 
= (-)5+i, and the fact that S = 0,1 are 
the only possibilities for a pair of spin 1/2 particles. 

What does this all mean for a model with an s-channel 
exchange particle coupling to Majorana bilinears? It 
means that among the basis fermion bilinears, the candi- 
dates are just the pseudoscalar ^i'fs'i' (0 ^), the scalar 
'I'^' jO++), and the axial vector '1'7^75\1> (1++). The vec- 
tor ^>jf^^/ (1 ), tensor ^'cr'*"* (2+~), and pseudotensor 
^'^75(7'"'^' (2^^) bilinears are C-odd and therefore dis- 
allowed. The only s-channcl particles which may couple 
to these candidate bilinears are the pseudoscalar, scalar, 
or axial vector. 

There is some subtlety associated with the s-channel 
exchange of an axial- vector. The axial- vector is an L = 1 
mode, and we have seen that this mode elicits a v"^ sup- 
pression in the rate. However, the exchange particle 
is off-shell (away form resonance) and so has a time- 
like pseudoscalar piece in addition to the axial three- 
vector piece. This pseudoscalar coupling is effectively 
dfi (^7''75^'). The weak interaction coupling of the pion 
to the axial vector current provides a familiar example 
of such a coupling. The axial current is not conserved, 
and so the pseudoscalar coupling is nonzero. One has 
9^(*7^75*) = 2im,*75«' - ^e^^-^^f^k^Kikykc^Xpik). 
The first term shows an m^-dependence in the ampli- 
tude, leading to (m^/Mx)'^ helicity-suppression of the 
L = piece, while the second term is the famous anoma- 
lous VVA coupling. It offers W^W~ and ZZ produc- 
tion (with momenta k, k and helicities A(fc), A(fc)), but 
at higher order aw = ,9v'/47r in the electroweak Wiyv 
or Zuu coupling gy- The linear combination of a v^- 



6 



suppressed L = 1 piece and a m^-suppressed L = piece 
to the rate from axial vector exchange was first noticed 
by Goldberg [13]. 

For s-channel exchange of a true pseudoscalar or scalar 
particle, there is no helicity suppression. In addition, it 
was shown several years ago [M] that helicity suppression 
in the axial-vector case may be avoided when the two- 
body final state is replaced by a three-body final state. 
In the work of [sj, a charged pair was produced and a 
photon was radiated from one of the charged particles. 
Radiation of the photon changes the Dirac structure of 
the current; it also allows one fermion to become virtual, 
with a large replacing its small mass. It was found that 
brehmsstrahlung from the final state was dominated by 
collinear and infrared emission, which left the emitting 
fermion nearly on-shell, and helicity suppressed. How- 
ever, photon emission from the accompanying f-channel 
particle did give an unsuppressed amplitude. In our 
work, it is a massive W or Z-boson that is radiated. 
There are no collinear of infrared singularites, and the 
virtual fermion is necessarily off-shell by ~ M^r■ It 
seems likely to us that this 2 — + 3 rate for the radiatively 
corrected s-channel axial- vector exchange will exceed the 
helicity-suppressed 2^2 rate by 37 ('^^)^, which is 
many orders of magnitude. We do not pursue this fea- 
ture of the axial-vector exchange further here, for the 
premise of the neutrino-only model, investigated in this 
work, is that the tree-level annihilation rate to neutrinos 
is unsuppressed. 

So far we have discussed s-channel exchange processes. 
We turn now to a brief discussion of f-channel exchange 
annihilation models. The implications of a Majorana 
XX PE^ir are best recognized by Fierz transforming the 
two fermion bilinears to "charge-retention" order, i.e., to 
a x-bilinear and a i^-bilinear. If the Fierz'd bilinears con- 
tain a pseudoscalar, there is no suppression of the rate. 
Otherwise, there is a rate suppression. If the Fierz'd 
bilinears contain an axial vector piece, it contributes to 
the rate a w^-suppresed piece, and a (mj^/M^)^ helicity- 
suppressed piece. (However, as just explained above, a 
W 01 Z radiated in the final state may well lead to a 
much enhanced rate.) 

It is illuminating to explain in this context the of- 
ten seen remark that LSP-annihilation has a helicity- 
suppressed rate to fermions. This is true for SUSY exten- 
sions of the SM. It is not true in general for models of dark 
matter. For SUSY extensions of the SM, the annihilation 



graphs consist of i-channel scalar exchanges, and from 
crossing the identical Majorana fermions, also u-channel 
scalar exchanges; in addition, there are scalar, pseu- 
doscalar and axial- vector s-channel exchanges. Fierzing 
the t- and u-channel scalar exchanges yields s-channel 
axial- vector bilinears [l8l | , with the concomitant helicity- 
suppressed L = contribution and w^-suppressed L — 1 
contribution to the annihilation rate. The only contribu- 
tions that are potentially large come from the s-channel 
pseudoscalars. However, the scalars and pseudoscalars 
are Higgs particles, whose Yukawa couplings gy to the 
SM fermion are all proportional to (ruf /vev), thereby 
giving the same effect as a true helicity suppression. Re- 
call that the Higgses are assigned the double burden of 
providing mass to the electroweak gauge bosons and to 
the fermions. The gy oc vev~^ relation can be traced 
back to mass generation of the gauge bosons, while the 
gy oc mf relation comes from mass generation for the 
fermions. There are other possibilities for the Yukawa 
coupling to neutrinos. To give one example, if the neu- 
trino mass is small as a result of a see-saw mechanism, 
then the coupling of the Higg doublet to the neutrino 
field will be ^ GeV/vev, similar to the coupling of the 
Higgs to most quarks. And importantly, a more general 
scalar or pseudoscalar field, not complicit in fermon mass 
generation, would couple with an arbitrary gy. 

The upshot of all this for our investigation is the follow- 
ing: With an s-channel scalar or pseudoscalar B-mcson 
exchange, there is no helicity suppression. The Yukawa 
coupling of the i?-meson to Dv (and to xx) is arbitrary. 
The scalar exchange proceeds in the L = 1 partial wave, 
which suppresses the xX annihilation rate by v'^. On 
the other hand, the pseudscalar exchange proceeds in the 
L = partial wave, with no suppression of the rate. 
These deductions from partial wave analysis hold true for 
Dirac or Majorana x- 

In the rate ratio that we investigate, any rate sup- 
pression factors out. Moreover, when our assumed scalar 
exchange is replaced with a pseudoscalar exchange, there 
is no change in the rate ratio. This is easily seen by noting 
that placement of an 175 in the amplitude of Eq. (7) just 
before the spinor v(ki) alters this amplitude by just the 
overall phase —i (neglecting the neutrino mass). Thus, 
our results as presented hold also for the unsuppressed 
L — s-channel exchange, all the way down in energy to 
the XX annihilation threshold. 
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